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ABSTRACT 


This  report  discusses  analytical  and  numerical  approaches  for  the  theoretical  analysis  of 
adhesively  tended  structures  having  various  tendline  quality  levels.  Since  the  adhesive 
strength  of  an  adhesively  tended  structure  is  the  critical  failure  problem  of  a  bond  under 
small  loads,  the  discussions  in  this  paper  concentrate  on  mathematical  models  which  use 
interfacial  layers  between  the  adherend  and  adhesive  layers  as  a  means  to  show  various 
degrees  of  tend  quality.  The  analytical  method  discussed  is  a  reflection  coefficient  calcula¬ 
tion  for  multi-layered  structures  in  the  frequency  domain.  The  numerical  approach  used  is 
the  finite  element  method  yielding  time  domain  results  for  adhesively  tended  structures  hav¬ 
ing  various  interfacial  layer  conditions.  Direct  comparisons  between  the  two  methods  are 
accomplished  in  the  frequency  domain,  the  time  domain  results  obtained  via  the  finite  ele¬ 
ment  method  being  appropriately  Fast  Fouier  Transformed  (FFT). 
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INTRODUCTION 


The  strength  of  an  adhesive  bond  depends  both  on  its  cohesive  and  adhesive  properties. 
Cohesion  is  primarily  a  function  of  the  type  of  adhesive,  its  mechanical  properties,  degree  of 
cure,  porosity,  and  thickness,  all  being  properties  of  the  adhesive  bond’s  volume.  Adhesion  is 
dependent  on  the  degree  of  molecular  attraction  and  interaction  between  the  adhesive  and 
adherends,  depending  primarily  on  surface  area  conditions  of  the  adhesive  bond.^ 

In  typical  metal  to  metal  adhesively  bonded  joints,  surfaces  are  cleaned  and  primed  prior 
to  bonding.  This  process  results  in  the  formation  of  an  interfacial  layer  typically  only  a  frac¬ 
tion  of  a  micrometer  thick. ^  Improper  surface  preparation  degrades  the  quality  of  this  layer 
and  ultimately  causes  adhesive  bond  failure. 

Researchers^'^’^’^®  have  tried  to  characterize  adhesive  bond  properties  using  ultrasonic  non¬ 
destructive  evaluation  time  and  frequency  domain  techniques.  In  most  of  the  research,  adhesively 
bonded  structures  were  modeled  as  three-layered  media;  i.e.,  as  adherend/adhesive/adherend 
models.  Experimental  measurements  are  made  in  the  time  domain,  transformed  into  the 
frequency  domain,  and  analyzed  for  changes  in  resonance  minima:  spacing,  width,  and  depres¬ 
sion  depth.  Others®’^  have  taken  into  account  the  interfacial  layer  and  typically  use 
obliquely  incident  ultrasonic  wave  methods  for  interrogation.  In  particular,  Rose^  has  clearly 
shown  that  obliquely  incident  shear  wave  propagation  aimed  toward  an  interface  can  produce 
a  significant  increase  in  sensitivity  in  comparing  good  versus  weak  interface  areas  in  an 
adhesively  bonded  structure. 

In  this  paper,  the  effects  of  changing  interfacial  layer  densities  are  examined  to  simulate 
various  adhesive  bond  quality  levels  using  normal  incident  longitudinal  waves.  Analytical  and 
numerical  methods,  specifically  finite  elements,  are  used  to  examine  the  ultrasonic  wave  reflec¬ 
tion  coefficients  and  mechanical  displacement  responses,  respectively,  of  an  adhesively  bonded 
structure  consisting  of  adherends,  interfacial  layers,  and  an  adhesive  layer  (five-layers).  Differ¬ 
ences  in  resonance  minima  location  and  spacing  are  correlated  with  interfacial  layer  parame¬ 
ters  and  suggest  ways  to  estimate  adhesive  bond  quality. 

ANALYTICAL  MODEL 


Background 

Brekhovskikh'^  showed  three  approaches  for  studying  acoustic  wave  propagation.  The  possi¬ 
ble  approaches  are:  (1)  by  solving  the  wave  equation  with  appropriate  boundary  conditions, 

(2)  by  means  of  an  input  impedance  approach,  analogous  to  transmission  line  theory,  and  (3) 
by  considering  the  multiple  reflections  that  occur  at  each  interface  of  the  modeled  structure. 
The  second  approach  is  used  in  this  paper  for  studying  an  adhesively  bonded  multi-layered 
structure. 

At  a  plane  boundary,  a  simple  interface  may  exist  between  two  different  materials,  having 
an  acoustic  plane  wave  incident  upon  and  propagating  from  the  upper  to  lower  material  with 
mechanical  impedances  of  Z2  and  zj  for  the  upper  and  lower  materials,  respectively.  The 
wave  reflection  coefficient  may  be  obtained  as  follows: 


zi  -Z2 

r  =  - 

Zi  -I-  Z2 


(1) 
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where  the  mechanical  impedance  of  any  material  is  defined  as  the  product  of  its  material  den¬ 
sity  and  material  acoustic  wave  velocity,  z  =  />  v. 

Brekhovskikh  shows  that  Equation  1  may  be  rewritten  in  the  form. 


Zin  ^2 
r  =  - 

^in  ^2 


(2) 


where  Zin  equals  zi  of  Equation  1  and  is  defined  as  the  impedance  of  the  reflecting  medium. 
This  leads  to  the  general  equations  for  calculating  the  input  impedance  and  reflection  coeffi¬ 
cient  that  can  then  be  used  successively  for  modeling  multi-layered  systems. 


The  geometry  for  plane  wave  propagation  in  multi-layered  systems  is  illustrated  in  Figure  1 
where  normal  acoustic  plane  incidence  is  assumed  here.  Both  the  top  and  bottom  media  are 
semi-infinite  in  length  with  an  acoustic  plane  wave  incident  on  the  top  boundary.  The  top 
media  is  numbered  n-fl,  the  bottom  media  is  number  1,  and  there  are  n-2  intermediate  layers 
numbered  2,3,4,... ,n. 


For  the  n^**  layer  the  input  impedance  is: 


zC")  = 

in 


Mm  n 


tan^n 


in 


Zn  -  ‘Z'  jn  tan^n 


(3) 


where  for  normal  incidence  =  2ji(/vn,  kn  is  the  wave  number,  dn  is  thickness, 

f  is  frequency,  and  Vn  is  the  wave  velocity. 

The  reflection  coefficient  for  the  n***  layer  is: 


r 


Zn+l 


+  Zn+1 


(4) 


which  has  the  same  form  as  Equation  1,  the  reflection  coefficient  for  a  one  boundary  situation. 

Model  Description 

A  typical  adhesively  bonded  structure  consists  of  two  adherends,  two  interfacial  layers  on 
each  side  of  the  adhesive  layer,  and  the  adhesive  layer  itself.  Aluminum  is  the  adherend 
material  for  this  case.  Generally,  aluminum  oxide  and  primer  layers  are  generated  in  the 
anodization  and  priming  processes  to  enhance  bond  quality.  Adherends  can  be  considered 
as  finite  thickness  layers,  or  when  short  ultrasonic  pulses  aie  used,  as  semi-infinite  spaces. 
Since  a  very  short  pulse  excitation  is  considered  in  the  application  presented,  comparisons 
may  be  justifiably  accomplished  between  both  these  adherend  situations.  The  aluminum  oxide 
and  primer  materials  are  modeled  as  one  interface  layer  both  above  and  below  the  adhesive 
material  due  to  the  very  small  thickness  dimensions  of  these  individual  components. 
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For  the  purposes  of  the  analytical  reflection  coefficient  model,  the  upper  and  lower 
media  are  considered  to  be  semi-infinite  in  length  while  the  interfacial  and  adhesive  layers 
are  finite  in  thickness.  The  adhesive  bond  and  interfacial  layers  are  shown  geometrically  in 
Figure  2.  There  are  five  regions;  three  are  intermediate  layers  d2,  d3,  and  d4  which  repre¬ 
sent  the  interfacial  layer  above  the  adhesive,  the  adhesive,  and  the  interfacial  layer  below  the 
adhesive,  respectively. 

The  input  impedance  for  this  five-layer  system  is  Zin(4),  obtained  through  successive 
applications  of  Equation  3,  being  recursively  substituted  in  Equation  4  yielding  the  reflection 
coefficient  of  the  layered  system  as  follows: 


r 


+  Z5 


(5) 


where  Z5  is  the  mechanical  impedance  of  the  5*'’  material. 

TTiis  reflection  coefficient.  Equation  5,  is  used  to  examine  the  effects  of  a  degrading  adhe¬ 
sive  bond  for  two  specific  cases.  A  perfect  bond  is  defined  with  the  interfacial  layer  having 
either  the  same  mechanical  impedance,  acoustic  wave  velocity,  and  material  density  as  the 
adherend  for  Case  1,  or  as  the  adhesive  for  Case  2.  Degrading  interfacial  layer  conditions 
are  modeled  by  varying  the  mechanical  impedance  of  the  interfacial  layers  through  decreases 
in  the  material  density  of  the  layer.  The  various  conditions  for  Cases  1  and  2  are  listed  in 
Tables  1  and  2,  respectively. 


NUMERICAL  MODEL;  FINITE  ELEMENT  APPROACH 


The  analytical  model  previously  described  is  valid  for  one-dimensional  wave  propagation 
analysis  and  for  situations  where  the  effects  of  the  second  and  third  dimensions  have  a  very 
small  consequence  on  the  mechanical  displacement  solution  of  the  acoustic  wave  equation. 
The  acoustic  wave  equation  for  damped  media  is  given  as  follows: 


dt2 


=  V. 


iijL  +  V  7, 

az2  *  ay2  ' 


dz^  dt 


+  V. 


u 

dy^dt 


+  F(t) 


(6) 


where  P  is  the  material  density  in  kg/m  ,  Vl  is  the  longitudinal  phase  velocity  in  m/sec,  Vj  is 
the  shear  phase  velocity  in  m/sec,  7l  is  the  mechanical  damping  parameter  in  N-s/m^  in  the 
z-direction,  7s  is  the  mechanical  damping  parameter  in  N-s/m^  in  the  y-direction,  and  F(t)  is 
the  applied  force  density  in  N/m^  on  the  structural  surface  as  a  function  of  time. 

The  finite  element  method  transforms  the  differential  equation  given  by  Equation  6  into 
an  integral  relatit  nship  :  nd,  finally,  into  an  algebraic  system  of  equations  which  may  be 
solved  on  a  digital  computer.  Details  of  this  procedure  will  not  be  given  here  but  may  be 
found  in  References  18  through  20.  A  two-dimensional  finite  element  analysis  is  given  on 
the  adhesively  bonded  structure  for  two  cases;  (1)  A  spatially  uniform  force  application  of 
finite  aperture  width  applied  by  a  transducer  (transducer  not  modeled)  into  an  adhesively 
bonded  structure  at  a  center  frequency  of  15  MHz  where  the  interfacial  layer  has  a  material 
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density  proportional  to  the  aluminum  adherend  (same  as  percent  of  bond  quality),  and  (2) 
same  as  Case  1  except  that  force  excitation  is  at  10  MHz  and  the  interfacial  layer  has  a  mate¬ 
rial  density  proportional  to  that  of  the  adhesive  material  (percent  of  bond  quality).  The  forc¬ 
ing  function  is  the  temporal  raised  cosine  function  as  follows: 

F  (t)  =  [1  -  cos  P  j  ]  cos  (2  ;r  fo  t)  0  <  t  s  ^  (7) 


where  fo  is  the  transducer’s  center  operating  frequency,  F(t)  is  the  time  domain  raised  cosine 
force  density  input  with  uniform  surface  application  to  the  structure. 

The  idea  of  using  the  Finite  element  method  is  that  it  solves  any  field  variable  of  interest 
at  discrete  nodal  locations  throughout  a  domain;  in  this  case  for  an  adhesively  bonded  struc¬ 
ture.  Mechanical  displacements  are  solved  in  this  case  using  this  method  for  the  five-layered 
structure  shown  in  Figure  3. 

The  choice  of  both  the  finite  element  spatial  discretization  level  and  the  chosen  time  step 
increment  for  simulating  the  mechanical  displacement  solution  must  be  done  appropriately  in 
order  to  ensure  valid  solution  computations.  For  the  purposes  of  the  spatial  discretization,  a 
spatial  sampling  of  8  nodal  points  per  center  frequency  wavelength  is  chosen.  By  the 
Nyquist  criterion,  a  minimum  of  two  samples  per  wavelength  is  required  to  represent  any 
signal.^^  Since  harmonics  of  the  fundamental  frequency  exist  in  the  mechanical  displacement 
response,  it  is  felt  that  Nyquist  spatial  sampling  up  to  the  fourth  harmonic  is  sufficient  for 
practical  purposes;  i.e.,  this  being  the  principal  reason  that  the  8-node  per  fundamental  fre¬ 
quency  wavelength  sampling  requirement  is  chosen. 

Since  the  time  domain  formulation  is  an  explicit  time  stepping  scheme,  the  time  step 
size  used  is  dependent  on  the  spatial  discretization;  in  this  case  the  8-node  per  wavelength  cri¬ 
terion.  The  matrix  equation  formulation  of  Equation  6  prior  to  time  discretization  is  given  by: 

[M]  {u} -f  [D]  {u} [K]  {u}  =  {F}  (8) 


where  [M]  is  called  the  mass  matrix,  [D]  is  the  mechanical  damping  matrix,  [K]  is  the 
mechanical  stiffness  matrix,  and  {u}  is  the  mechanical  displacement  vector  with  its  associ¬ 
ated  derivatives.**’^® 

The  finite  difference  time  discretization  of  the  second  and  first  derivatives  of  mechanical 
displacement  is  accomplished  by  Equations  9  and  10: 

Central  Difference  Approximation: 

(u,  +  At-2«,  +“,_Atl  (9) 


where  U{  is  the  central  difference  approximation  for  the  second  time  derivative  of  u,  t  is  the 
time  step  size,  and  u,.^At’  "t*  nodal  values  of  mechanical  displacements  at  the 

three  consecutive  time  instants,  t-t-At,  t,  and  t-At. 
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Backward  Difference  Approximation: 


"t  -  At 


(10) 


where  Ut  approximates  the  first  time  derivative  of  u  at  time  t. 

The  final  algorithm,  given  by  Equation  11,  used  for  the  mechanical  displacement  calcula¬ 
tions  combining  finite  elements  in  the  spatial  domain  and  finite  differences  in  the  time 
domain  is  given  by: 


At 


“t  +  A.  -  ‘  t  -  “t)  +  [D](u,-u,_^,)}  -I-  2u,  - 


(11) 


1 A 

where  Mjj  is  the  diagonalized  mass  matrix. 

Using  Fourier  Stability  Analysis,  the  stability  relationships  between  the  spatial  and  tempo¬ 
ral  discretization  levels  may  be  derived.  For  the  undamped  acoustic  wave  propagation  prob¬ 
lem,  the  final  stability  relationship  for  the  undamped  case  is: 

Ay  Az 

At<  -  .  ■  (12) 

^max  y  ^  y  +  A  z^ 


where  At  is  the  time  step,  Ay  is  the  spatial  discretization  in  the  y-direction,  Az  is  the 
spatial  discretization  in  the  z-direction,  and  Vmax  Ihe  maximum  phase  velocity  in  the 
structure. 

Tables  3  and  4  give  all  of  the  information  on  the  resulting  mesh  discretization  levels 
including  the  required  spatial  and  time  step  requirements  for  stable  algorithms  for  Cases  (1) 
and  (2)  previously  described. 

For  Case  1,  where  a  transducer  center  frequency  of  15  MHz  was  used,  and  the  adhesive 
interface  was  initially  modeled  like  aluminum,  the  required  time  step  size  was  1.060  nanosec¬ 
onds  and  the  actual  time  step  used  was  1.000  nanosecond. 

For  Case  2,  where  a  transducer  center  frequency  of  10  MHz  was  used,  and  the  adhesive 
interface  was  initially  modeled  like  the  adhesive,  the  required  time  step  size  was  3.280  nano¬ 
seconds  and  the  actual  time  step  used  was  1.000  nanosecond. 

Since  the  adhesively  bonded  structure  is  rectangular  in  shape  and  the  force  excitation 
is  applied  along  the  center  of  the  top  layer,  the  condition  of  symmetry  may  be  advanta¬ 
geously  applied.  Only  one-half  of  the  structure  needs  to  be  modeled,  either  the  right  or 
left  half,  by  assuming  a  zero  displacement  condition  in  the  y-direction  along  the  symmetry 
axis  (z-direction).  This  condition  is  shown  in  Figure  4  using  the  right  half  version  of  the 
model. 


5 


NUMERICAL  AND  ANALYTICAL  RESULTS  DISCUSSION 


Numerical  results  calculated  via  the  finite  element  method  are  time  domain  results  which 
are  found  at  every  discrete  time  step  for  the  time  interval  simulation  desired  at  all  the  nodal 
locations  in  the  given  domain  modeled;  in  this  case  the  adhesively  bonded  structure.  One  ob¬ 
servation  node  is  chosen  at  the  top  center  node  of  the  physical  structure,  or  equivalently,  the 
upper  left  node  of  the  symmetric  model  for  the  finite  element  simulation  (Figure  4).  The 
time  and  frequency  domain  results  presented  are  based  on  this  particular  nodal  location. 

The  comparisons  between  the  analytical  and  numerical  results  are  done  totally  in  the  fre¬ 
quency  domain.  All  time  domain  results  obtained  for  the  observation  node  are  shown  for  the 
time  interval  of  the  first  echo  pulse  only  since  this  timeframe  gives  the  pertinent  information 
about  the  bondline.  These  results  are  then  Fast  Fourier  Transformed  (FFT)  into  the  fre¬ 
quency  domain.  The  overall  frequency  range  is  studied  for  the  case  of  a  100%  quality  level 
bond  for  Cases  1  and  2.  For  Case  1,  the  effective  bondline  thickness  is  70  microns,  and  a 
resonant  frequency  (amplitude  minimum)  at  15  MHz  is  detected  analytically.  For  Case  2,  the 
effective  bondline  thickness  is  84  microns,  and  a  resonant  frequency  (amplitude  minimum)  at 
12.5  MHz  is  detected  analytically.  These  minima  correspond  exactly  to  the  one-dimensional 
analyses  for  bondlines  with  such  thickness  dimensions.  Both  the  reflection  coefficient  results 
and  the  displacement  output  frequency  responses  show  these  distinct  minima  very  clearly  in 
Figures  5  and  7  (for  the  100%  quality  level  simulation)  for  Case  1,  and  Figures  10  and  11 
(for  the  100%  quality  level  simulation)  for  Case  2. 

The  real  advantage  of  using  a  numerical  method  is  in  its  ability  to  solve  problems  where 
analytical  methods  are  not  possible.  It  is  always  necessary,  however,  to  compare  with  prob¬ 
lems  which  can  be  solved  analytically  in  order  to  gain  confidence  in  the  chosen  numerical 
model.  It  is  for  this  reason  that  the  case  of  a  good  bond  (a  100%  quality  bond)  excited  with 
a  normal  incidence  pulse  is  performed  and  its  results  are  shown  initially.  FFTs  are  performed 
on  the  numerical  time  domain  results  shown  in  Figures  8  and  12  and  the  corresponding  fre¬ 
quency  domain  results  of  the  analytical  and  numerical  results  are  appropriately  overlayed  in 
Figure  14.  There  is  an  excellent  comparison  for  Cases  1  and  2  between  the  numerical  and 
analytical  analyses,  especially  in  the  vicinity  of  frequency  minima  locations.  The  comparisons 
for  the  overall  frequency  range,  however,  are  not  one-to-one  although  they  are  good.  This  is 
attributed  to  several  reasons.  The  finite  element  method  employed  for  the  numerical  approxi¬ 
mation  here  is  a  two-dimensional  versus  a  one-dimensional  model.  Additionally,  mode  cou¬ 
pling  exists  between  the  longitudinal  and  shear  responses  of  the  adhesively  bonded  structure. 
This  is  true  although  the  present  structure  is  driven  with  a  force  distribution  only  in  the  longi¬ 
tudinal  direction.  This  mode  coupling  is  accounted  for  in  the  different  magnitude  values  of 
the  frequency  spectra  of  the  output  responses  of  the  analytical  versus  numerical  results.  It  is 
noteworthy  to  observe  that  the  frequency  values  of  the  numerical  results  are  consistently 
higher  than  the  analytical  ones.  This  is  characteristic  of  a  finite  element  formulation  using  an 
explicit  time  stepping  scheme  based  on  the  central  difference  operator.*^ 

The  adhesively  bonded  structures  with  partially  good  quality  bondlines  are  investigated  for 
10%,  20%,  and  50%  bond  quality  levels.  Tlie  finite  element  time  domain  results  are  shown 
in  Figures  8  and  12,  transformed  into  the  frequency  domain  in  Figures  9  and  13,  and  compari¬ 
sons  are  made  for  the  minima  frequency  results  between  the  analytical  reflection  coefficient 
(Figure  6)  and  numerical  results  for  Cases  1  and  2.  These  minima  results  are  shown  for  the 
different  levels  of  bond  quality  in  Table  5.  Note  that  the  error  between  the  frequency  min¬ 
ima  is  always  less  than  about  2%  demonstrating  that  the  numerical  approach  is  certainly  a 
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credible  approach  for  analyzing  adhesively  bonded  structures  with  bondlines  of  variable  quality 
levels.  Additionally,  a  correlation  may  be  reached  by  inspecting  the  widths  of  frequency  spec¬ 
tra  in  the  vicinity  of  the  amplitude  minima  locations.  These  widths  or  "bandwidths"  appear  to 
be  somewhat  proportional  in  size  to  the  quality  of  the  adhesive  bondline  being  analyzed;  i.e., 
being  widest  for  the  higher  quality  bonds.  This  was  consistently  found  not  only  in  the  compar¬ 
isons  between  numerical  and  analytical  results  presented  here  but  also  in  a  theoretical  and 
experimental  development  made  by  Chang.^  This  showed  that  spectral  depression  bandwidths 
may  be  indicators  of  bond  strength  as  well  as  being  related  to  the  acoustic  interfaces  of  a 
bonded  joint.^ 


CONCLUSIONS 

This  paper  demonstrates  that  the  finite  element  method  is  a  feasible  means  of  evaluating 
adhesively-bonded  structures  where  the  quality  of  the  bondline  enters  directly  into  the  analy¬ 
sis.  Since  adhesive  failures  occur  with  a  higher  probability  than  cohesive  failures,  an  interfa¬ 
cial  layer  was  introduced  to  simulate  weak  bondline  conditions  in  the  modeled  structure.  The 
comparison  of  numerical  versus  analytical  methods  shows  that  the  detection  of  frequency  min¬ 
ima  by  numerical  methods  for  the  five-layered  model  simulated  for  various  bond  quality  levels 
indeed  gives  very  similar  resonant  frequency  results  to  analytical  methods.  Additionally,  the 
usage  of  interfacial  layers  as  a  means  to  simulate  weak  conditions  at  a  bondline  surface  area 
junction  appears  to  be  an  excellent  way  to  model  the  adhesive  failure  conditions  usually 
encountered.  Characteristics  which  occur  in  the  frequency  domain  for  various  quality  bondline 
conditions  verified  by  the  finite  element  model  are:  (1)  a  reduction  in  the  resonant  fre¬ 
quency  location  for  reduced  bondline  quality  level,  (2)  a  reduction  in  the  bandwidth  around 
the  frequency  minima  location  for  reduced  bondline  quality  level,  (3)  excellent  comparisons 
with  analytical  reflection  coefficient  -  lodel  (two-dimensional  versus  one-dimensional  analyses), 
and  (4)  verification  of  the  interfacial  layer  since  it  is  qualitatively  consistent  with  the  theoreti¬ 
cal  and  experimental  discussion  given  by  Chang.^^ 


Table  1 .  CASE  1  -  ANALYTICAL  PARAMETER  LIST 


Velocity  (v) 

(m/s) 

Density  (p) 
(kg/m^ 

Impedance  (z) 
(xIO®  kg/m^  sec) 

Thickness  (d) 

(m) 

Adherend 

6370 

2710 

17.25 

Semi-infinite 

Adhesive 

2100 

1120 

2.35 

70.0  X  10"® 

Interface  Layers:* 

Model-a  100% 

6370 

2710 

17.25 

7.0x  lO'® 

Model-b  50% 

6370 

1360 

8.66 

7.0  X  10"® 

Model-c20% 

6370 

540 

3.44 

7.0  X  10"® 

Model-d  10% 

6370 

270 

1.72 

7.0x10^ 

*l-ercentage  of  adherend  density  corresponding  impedance  calculated  using  z=pv 
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Table  2.  CASE  2  -  ANALYTICAL  PARAMETER  LIST 


Velocity  (v) 

(m/s) 

Density  ip) 
(kg/m^ 

Impedance  (z) 
(xIO®  kg/m^  sec) 

Thickness  (d) 

(m) 

Adherend 

6370 

2710 

17.25 

Semi-infinite 

Adhesive 

2100 

1120 

2.35 

70.0  X  10"® 

Interface  Layers;* 

Modei-a  100% 

2100 

1120 

2.35 

7.0x10'® 

Modet-b  50% 

2100 

560 

1.18 

7.0x10'® 

Modei-c20% 

2100 

220 

0.47 

7.0x10'® 

Model-d  10% 

2100 

110 

0.23 

7.0x10'® 

^Percentage  of  adhesive  density,  corresponding  impedance  calculated  using  z=pv 


Table  3.  CASE  1  -  NUMERICAL  PARAMETER  UST 


Thickness 

(m) 

Vl 

(m/s) 

Vs 

(m/s) 

#Z  Elements 

Al/8 

(m) 

Is/B 

(n.) 

At 

(nsec) 

AZ 

(m) 

AY/AZ 

3.17  x  10'® 

6370 

3110 

123 

5.3  X  10'® 

2.59  X  10® 

2.861 

2.58x10® 

0.998 

7.00x10^ 

6370 

3110 

1 

5.3  X  10® 

2.59x10® 

1.060 

7.00x10'® 

3.68 

7.00  X  10'® 

2100 

1050 

9 

1.75x10® 

8.75  X  10^ 

3.542 

7.77x10'® 

3.31 

7.00x10'® 

6370 

3110 

1 

5.3  X  10® 

2.59x10® 

1.060 

7.00  X  10^ 

3.68 

4.77  X  10'® 

6370 

3110 

185 

5.3  X  10® 

2.59x10® 

2.861 

2.58x10® 

0.998 

#DOF  =  316160 

nz  =  320 

ny  = 

494  At = 

1.060  nsec 

Table  4.  CASE  2  -  NUMERICAL  PARAMETER  UST 

Thickness 

(m) 

Vl 

(m/s) 

Vs 

(m/s) 

#Z  Elements 

Al/8 

(m) 

As/8 

(m) 

Ar 

(nsec) 

AZ 

(m) 

AY/AZ 

3.17  X  10'® 

6370 

3110 

82 

7.96  X  10  ® 

3.88  X  10® 

4.296 

3.87x10® 

1.000 

7.00x10'® 

2100 

1050 

1 

2.625  X  10  ® 

1.313x10® 

3.280 

7.00x10"® 

5.531 

7.00x10® 

2100 

1050 

6 

2.625  X  10® 

1.313  X  10"® 

5.316 

1.166x10"® 

3.319 

7.00x10'® 

2100 

1050 

1 

2.625  X  10  ® 

1.313x10® 

3.28C 

7.00x10"® 

5.531 

4.77x10"® 

6370 

3110 

123 

7.96  X  10® 

3.88x10® 

4.304 

3.882x10® 

09974 

#DOr  =  140812 

nz  =  241 

ny  = 

329  At 

=  3.280  nsec 
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Table  5.  FREQUENCY  MINIMA  COMPARISON 


Frequency  Minima  Locations 

Case/Model 

Density  of 

Interface  Layers 
(kg/m^ 

Analytical 

(MHz) 

Numerical 

(MHz) 

Case  1 

Model-a 

2710 

15.00 

14.93 

Modei-b 

1360 

14.80 

14.73 

Model-c 

540 

14.35 

14.31 

Model-d 

270 

13.75 

13.72 

Case  2 

Model-a 

1120 

12.50 

12.40 

Model-b 

560 

10.80 

10.80 

Modet-c 

220 

8.20 

8.12 

Model-d 

110 

6.20 

6.16 

I 

I  TRANSMITTED 


Figure  1.  Multi-layered  structure  model  used  for  calculation  of  reflection  coefficient. 
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Figure  2.  Analytical  representation  of  adhesively  bonded  structure 
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Figure  3.  Adhesively  bonded  structure  geometry  with  finite  aperture 
transducer  pulsed  excitation  (adherends  are  aluminum). 
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Figure  4.  Possible  finite  element  mesh  of  five-layered  structure  exploiting  mechanical 
symmetry,  shear  displacements  set  to  zero-value  on  symmetry  axis. 


case  1.  1007.  analytical 


FREQUENCY  (Wz) 


Figure  5  Analytical  reflection  coefficient  Case  1  (good  bond) 
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RELATIVE  AMPLITUDE 


case  1.  analytical 


FREQUENCY  Ofb) 

Figure  6.  Analytical  reflection  coefRcients  for  Case  1 .  Bond  quality 
levels:  10%,  20%,  50%,  100%  •  13.5  s  f  s  15,5  MHz. 
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(a)  TRANSOJCER  SPECTRA 


H 

1 

m 

f 

/ 

y  > 

\ 

•  S  U  IS  ai  25  « 

HOBOr  Mkt 


<b)  NBUETICN  SCCIRA 


(c)  OUTPUT  SPECTRA 


FREQUENCY  (mz) 


Figure  7  Analytical  frequency  domain  analysis  for  good  bond 
(100%),  Case  1 :  (a)  Fourier  transform  of  raised  cosine  with  fo  =  15 
MHz,  (b)  Reflection  coefficient  of  adhesively  bonded  model,  and 
(c)  Relative  amplitude  of  transformed  mechanical  displacement 
response. 
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cae  1.  nuKrical 


Figure  8  Time  domain  mechanical  displacement  responses,  Case  1 : 
(a)  100%  txjnd  quality  level,  (b)  50%  bond  quality  level,  (c)  20%  bond 
quality  level,  and  (d)  10%  bond  quality  level 
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(ELAnVE  MTUIUDE 


RELATIVE  AMPLITUDE 


case  2.  analytical 


FREQUENCY  (MHz) 


Figure  10  Analytical  reflection  coefficients  Case  2,  5  <  f  s  13  MHz 
for  10%,  20%,  50%,  and  100%  bond  quality  levels  • 
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sPEnm 


CIS*  2.  nuaerical 


Figure  12  Time  domain  mechanical  displacement  responses,  Case  2 
(a)  100%  bond  quality  level,  (b)  50%  bond  quality  level,  (c)  20%  bond 
quality  level,  and  (d)  10%  bond  quality  level. 
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FREQUENCY  (MHz) 


Figure  13  Frequency  domain  mechanical  displacement  responses, 
Case  2,  for  10%,  20%,  50%.  and  100%  bond  quality  levels 


(a)  case  1.  analytical /ncunerical 
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(b)  case  2.  analytical/numerical 
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Figure  1 4  Frequency  domain  comparisons  of  analytical  versus  numeri¬ 
cal  approaches  for  good  bond  (100%  bond  quality  level):  (a)  Case  1, 
interfacial  layer  similar  to  adherend,  and  (b)  Case  2,  interfacial  layer 
similar  to  adhesive 
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